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■ Abstract 

We study correlation functions of elementary fermions in strongly interacting field 
theories using the AdS/CFT correspondence. This correspondence generically 
associates bulk fields to composite operators in field theory. We modify the holo- 
graphic prescription in order to obtain correlators that correspond to fermonic 
single-particle excitations by introducing a dynamical fermionic source localized on 
a UV brane in a holographic background. We work out the conditions when these 
correlators obey the zeroth frequency sum-rule satisfied by angle-resolved photo- 
emission spectroscopy (ARPES) and are thus directly relevant to the AdS/CMT 
correspondence. To illustrate our techniques, we study field theories at zero chem- 
ical potential with an arbitrary dynamical exponent z, i.e., the Lifshitz invariant 
conformal field theories, including the usual relativistic case z = 1. 
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1 Introduction 



The AdS/CMT program, see [T] and references therein, is a useful tool for the 
description of strongly interacting condensed-matter systems that grew out of the 
AdS/CFT correspondence P] recently. In its first instance it entails a duality 
between a gravity theory in a {d + l)-dimensional anti-de Sitter (AdS) bulk space- 
time and a strongly interacting conformal quantum field theory (CFT) living on 
the d- dimensional boundary of the AdS space. Perturbations away from the exact 
quantum critical point described by the conformal field theory can be realized, for 
instance, by considering a black brane in the anti-de Sitter space-time leading to a 
nonzero temperature for the theory on the boundary. Furthermore, a black brane 
charged with respect to an electro-magnetic gauge field corresponds to adding a 
chemical potential to the CFT. Studies of fermionic correlation functions in these 
systems lead to interesting Fermi and non- Fermi- like behavior [31 HJ El E] . 



T 



QCP 



Figure 1: The phase diagram of a particle- hole symmetric semi- metal. The tem- 
perature of the semi-metal is denoted by T and the chemical potential, which de- 
termines the nature and amount of doping of the semi-metal, by fj,. The dashed 
lines denote the smooth crossover between the classical and quantum behavior of 
the normal phase of the semi-metal. The latter is governed by the quantum critical 
point (QCP) at zero temperature and zero chemical potential. 



To understand that the AdS / CFT correspondence may also be useful for con- 
densed-matter physics, let us revisit the phase diagram of a particle-hole symmet- 
ric semi-metal. It is shown in figure 1, in a manner convenient for this discussion. 
Plotting the phase diagram in this way clearly shows that the behavior of the 
semi-metal is governed by a quantum critical point at zero temperature and zero 
chemical potential. Within the context of the AdS/CFT correspondence, the con- 
formal field theory associated with this quantum critical point has a particular 
important feature, namely that under a scale transformation the time direction 
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scales equally fast as the spatial directions. This is due to the relativistic invariant 
nature of the anti-de Sitter background, and from the point of view of condensed- 
matter physics it corresponds to a quantum critical point that has a dynamical 
exponent z = 1. Interestingly, graphene presents a concrete realization of this 
situation. However, in condensed-matter physics one more often deals with the 
situation that under a scale transformation the time direction scales twice as fast 
as the spatial directions, i.e., the dynamical exponent is z = 2. Examples are 
double-layer graphene or gapless semi-conductors. The latter are usually not ex- 
actly particle-hole symmetric, but in a first approximation this will not change the 
physics qualitatively. A system with z = 2 that breaks particle-hole symmetry in 
a more fundamental manner is given by an atomic Fermi mixture at unitarity. 

As mentioned above, anti-de Sitter space-time does not obey this so-called 
Lifshitz scaling but instead a relativistic scaling with z = 1. Therefore, the usual 
AdS/CFT correspondence has to be modified to be useful for this purpose. One 
possibility is to geometrically realize the full Schrodinger symmetry of the quantum 
critical point as in [71 E] which, however, has been shown in P [ID] to lead to 
different thermodynamical properties than fermions at unitarity. Another way 
to study the duality with an arbitrary dynamical exponent z is to consider a 
simple generalization of the anti-de Sitter space-time that exhibits the appropriate 
dynamical scaling and is now known as the Lifshitz space-time [TT]. Here, particle- 
hole symmetry is still an exact property of the conformal field theory on the 
boundary, just like in the anti-de Sitter case. This is the background that we use 
in this worl(0. 

1.1 Single-particle correlators and ARPES sum-rules 

From the condensed-matter point of view we are interested in calculating the 
single-particle correlation function, or Green's function, which for electronic sys- 
tems is directly observable by angle-resolved photo-emission spectroscopy (ARPES), 
or for ultracold atoms by radio-frequency (RF) spectroscopy. In particular, the 
single-particle Green's function G{k, u) differs from the retarded correlation func- 
tion of a composite operator in that the former obeys a special sum-rule which is a 
direct consequence of the canonical (anti) commutation relations of the elementary 
fields. It reads 



^ These backgrounds also exhibit subtleties as they are not geodesically complete and give 
rise to diverging tidal forces at the origin of the spacetime, see [12] for a very recent discussion. 
We will avoid these issues in this work by regulating the theory by turning on an infinitesimal 
black-hole horizon around the singularity. 
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A bulk fermion field \1/ in the AdS / CFT correspondence would typically correspond 
to composite operators, schematically of the form Tr (0 ■ • • (pipcj) ■ ■ ■ 0) where and 
ip represent bosonic and fermionic elementary fields, respectively, in the adjoint 
representation of the dual gauge theory. These operators satisfy modified sum- 
rules that differ from ( II .ip . see [TB] for a recent discussion. Intuitively it is clear 
that f ll.ip will not be obeyed because of the UV divergence that arises from placing 
more than one operator at the same point in space-time. 

To see precisely how the elementary sum-rule fails with the standard prescrip- 
tion of the AdS/CFT correspondence, it is important to realize that the Dirac 
equation in the [d + l)-dimensional bulk contains an undetermined mass param- 
eter M that is related to the scaling dimension of the dual fermionic operator. If 
we consider the case of d even, then we can split any fermionic operator into its 
chiral components 0±. The scaling dimensions for these operators are then given 
by 

d + z — 1 , , 

At = ^ TM. (1.2) 

From the AdS / CFT correspondence we obtain, as shown explicitly in this paper, 
that the Green's function of the operator O- satisfies Go_ (0, oj) oc uP'^l^ at zero 
temperature and at momentum /c = 0. With this form of the Green's function, the 
only way to satisfy the requirement f ll.ip in the usual AdS/CFT correspondence 
is to choose the mass of the bulk fermion dual to the operator O- as M = —z/2, 
in the appropriate units set by the radius of the Lifshitz background. Naively 
this would correspond to the trivial case of a free theory. However, a detailed 
analysis, which is somewhat more difficult in the special case M = —z/2, shows 
that there are logarithmic corrections that yield a Green's function that obeys the 
exact Dyson equation 

G{0,u) = (1.3) 

oo — L(U, OJ) 

with a self-energy of the form S(0, cj) ~ cjlogcu. This situation indeed satisfies 
the zeroth frequency sum rule, but instead it spoils the causality of the theory 
that is signaled by a violation of the Kramers-Kronig relation. Put differently, the 
logarithmic corrections to the self-energy lead to the appearance of a (Landau) pole 
in the upper half of the complex frequency plane, which is physically unacceptable. 

It therefore appears that the standard prescription of the AdS/CFT corre- 
spondence does not allow for the calculation of the single-particle Green's function 
that is of crucial importance in condensed-matter physics as well as for strongly 
correlated systems. As mentioned above, this is of course expected, as the stan- 
dard AdS/CFT correspondence calculates Green's function of composite operators, 
not of the single-particle excitations that satisfy canonical ( ant i-) commutation re- 
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lation^. In view of this situation, the main aim of this paper is to point out 
a modified prescription that does allow for the extraction of a single-particle 
Green's function that obeys both the zeroth frequency sum rule and the Kramers- 
Kronig relation. We will find that this prescription can in fact be applied when 
-z/2 < M < z/2. 

The key ingredient of our construction is to introduce a UV cut-off surface close 
to the boundary, on which the dynamics of the elementary fermion interacting with 
the CFT is defined. This introduces a UV scale such that, at zero temperature 
and at zero spatial momenta, there is still a dimensionful scale at the critical point 
that allows for non-trivial single-particle Green's functions. The idea is to apply 
the alternative quantization where one integrates over the boundary value of the 
bulk field, on the UV cut-off surface in the presence of a kinetic term for this bulk 
value, instead of fixing it to be a stationary source. This procedure computes the 
two-point function of the dynamical source, which we propose to correspond to 
the elementary fermion in the corresponding field theory. 

As we will show, the new prescription leads to the following retarded single- 
particle Green's function for a fermion coupled to the CFT at vanishing spatial 
momenta, temperature and chemical potential 

G'M = — ~i , . ■ (1-4) 

Here, ^ is a dimensionful coupling constant related to the position of the UV-brane. 
We will show that with this prescription both ARPES sum rules and Kramers- 
Kronig relations are satisfied. 

The aforementioned restriction —zjl < M < z/2 can easily be understood 
from the field theory point of view. In view of f ll.2p we observe that M = —z/2 
corresponds to the unitarity bound for the operator (9_ of dimension A+. Below 
this value, the expression (11. 4p has singularities in the upper half plane, thus 
both the Kramers-Kronig relations and the sum-rule fails. On the other hand, 
recalling that Go_iuj) oc w^^/^, one finds that, when M > z/2 the interaction 
term becomes irrelevant in the IR, but instead it modifies the UV. Therefore 
the sum-rule hence the unitarity fails (although one can satisfy Kramers-Kronig 
relations with an appropriate sign of g in this case). This restricts the range of the 
operator coupled to the elementary field to lie between —z/2 < M < z/2. We find 
the corresponding restrictions on the gravitational side. As usual, the lower bound 
M > —z/2 corresponds to the requirement of finite energy for the bulk fluctuation 
corresponding to the operator O. The upper bound M < z/2 on the other hand. 



^As an example of such an elementary field, one can consider the J\f — -i super Yang-Mills 
theory with gauge group U{N) and think of the trace of one of the four canonical Dirac fermions 
in the theory. In the SU{N) theory this would be trivial, whereas it yields an elementary 
gauge- invariant fermionic operator in the U{N) theory. 
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corresponds to requiring that the bulk field is normalizable near the boundary. If 
M > z/2, one cannot apply the aforementioned alternative quantization, instead 
one has to use the regular quantization, which would yield the inverse of fll.4p . 

The introduction of extra elementary degrees of freedom in this manner is not 
novel: Contino and Pomarol [H] considered the possibility of having a kinetic 
term for a fermion on a cut-off surface, hence a dynamical source field, and have 
analyzed its consequences for the dual gauge theory in the context of beyond- 
the-standard-model physics. However, our prescription in section 12.21 differs from 
[Ti] in that, although we stay within the unitarity bound, we are able to generate 
a Green's function of the form fll.4p by considering the Green's function of the 
dynamical source instead of the operator (9_ that couples to it. Perhaps the most 
recent and relevant work that has considered introduction of elementary fermions 
in the context of AdS/CFT is [16], building upon the ideas in [15]. In this "semi- 
holographic" method, an additional elementary fermion field is introduced and 
coupled to an operator (9_ in the CFT. The Green's function of the elementary 
field is obtained by field-theory methods, in terms of the correlator of which 
results in an expression similar to (11.41) . In a sense, the method that we advocate 
in this paper can be viewed as a derivation of the semi-holographic method. We 
propose a derivation of the expression in (11. 4p by holographic methods, with no 
ad hoc introduction of an extra degree of freedom. In this paper, we focus here on 
the more universal sector with zero chemical potential and leave the generalization 
to nonzero chemical potential to future work. 

The remainder of this work is organized as follows. In section 2 we study a 
Lifshitz background in the presence of a neutral black-brane and derive the cor- 
responding Dirac equation. Furthermore, we present our prescription to compute 
the Green's function for elementary fermions, for all values of the critical exponent 
z. In section 3, we present and analyze the results of the above-mentioned pre- 
scription. Note that since we are considering a neutral black-brane, we describe 
the normal, but quantum critical, state of the system at zero chemical potential 
where particle- hole symmetry is exact. Finally, we end in section 4 with some 
conclusions and an outlook for future work. 



2 Generalities 

2.1 Lifshitz backgrounds 

We start by considering the action for a gravity system coupled to a gauge field 
with field strength Ff^,y and a scalar field (p in d + 1 space-time dimensions. 



16nGd+i 



R-A-\{d,<f>W<P)-\e"fF,,F^'' 



(2.1) 
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The conventions for the space-time signature is ( — h • • • +), and the sign of the 
cosmological constant A is negative for negatively curved space-times. Moreover, 
we use units in which h = c = 1 and Gd+i is the dimensionless Newton couphng 
constant. All length scales are measured in terms of the AdS or Lifshitz radius. 
A black-brane solution that asymptotes to Lifshitz spacetime is given by \T7] 

ds' = - V\r)r'^dt' + r'dx' , V\r) = I - (-) • (2.2) 

The radial coordinate r runs from the boundary at r = oo to the horizon of the 
black brane at r = r/j, and the temperature of the black brane is obtained by 
demanding the absence of a conical singularity at r/j which leads to 

T = ^ ^ ir,y . (2.3) 

Note that the units of temperature are such that Boltzmann's constant /cb = 1 
and that the metric (12. 2p enjoys the following Lifshitz isometry 

r^/r, t^l~H, x^l-^x, T ^ FT . (2.4) 

For the field strength F^,^ and the scalar field 0, the solution to the equation of 
motion is given in terms of an arbitrary constant / 

Frt = fr'^"-^ , e^^ = 2{z-l){z + d-l) . (2.5) 

Furthermore, the dilaton exponent A and cosmological constant A are fixed to be 

A = -y^^^, A=-(^ + ci-l)(z + ci-2). (2.6) 

Notice that in the limit z — )■ 1, the scalar field decouples. The theory then reduces 
to the neutral black-brane solution in anti-de Sitter space-time with isotropic scal- 
ing properties [T71 HHj. Hence our results also allow for the case z = 1. For z > 1, 
note that the dilaton diverges at the boundary, and also the field strength diverges 
in such a way that a charge and chemical potential cannot be defined. To introduce 
a chemical potential, another gauge field needs to be considered whose solution has 
proper boundary conditions [IB]. For the purpose of the present paper, namely a 
study of spectral-density functions at zero chemical potential, this is however not 
needed. 

The diverging asymptotic behavior is thought to be a problem for holography 
and a proper treatment of the holographic renormalization procedure is therefore 
needed. This has been initiated in [2U | [2T|. though not for Lifshitz space-times 
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supported by a scalar field like in the present paper. However, this problem is not 
severe for the applications we have in mind for the following reason. The fermions 
that we introduce below do not couple to the dilaton and are not charged under 
the gauge field. They only couple to the metric, which has a good asymptotic 
behavior. Therefore, the fermionic boundary Green's functions we obtain in this 
paper are well defined, as we show explicitly. 

2.2 Prescription to calculate Green's functions for elemen- 
tary fermions 

In order to study fermions in the dual boundary theory, we add the following 
action to the gravity system given by (12. ip . 

Sj[^] = igf j d'^^^x^/^ Q^^^ - ^^^^ - M^^^ + Sa[^] . (2.7) 

Our notation here is as usual: = ^^T^, f = V^Ca^V^ with V^, = d^ + \uj^,gh^^, 
and F- are the {d+ l)-dimensional gamma matrices, where F- is antihermitean in 
our conventions. Furthermore, are the vielbeins and w^ab is the spin con- 
nection, for which indices are raised and lowered with the fiat metric r]ab = 
diag (— 1, +1, . . . , +1). Finally, we defined F— = i[F-, F-]. Next, since we are 
interested in a space-time with boundaries, we introduce an action Sg^if] to make 
the variational problem well-posed [22| [H] . More concretely, we consider bound- 
aries at the horizon r = rt a.s well as at a cut-off surface located at r = tq. But, as 
we will see below, due to the fact that y/g^ vanishes at the horizon we only need 
to add a boundary term to the action for the cut-off surface at r = rg. 

Next, let us define components \E'± of the Dirac fermion \l/ using the gamma 
matrix F- as follows 

vi/±^l(l±F^)vl/, F^vl>^ = ±v^±, (2.8) 

where for d being odd F- is the radial gamma matrix, while for d being even it 
is the boundary chirality operator. In both cases, we choose it to be hermitian. 
Note that a Dirichlet boundary condition for either \Ef_|_ or \E'_ can be imposed 
consistently, however, one cannot impose both 5\E'+ = and 5\E'_ = at the same 
time because the Dirac equation is of first order in derivatives. Mixed boundary 
conditions are also possible, and for the case when Lorentz symmetry is broken on 
the boundary see [23]. Here, we focus on pure Dirichlet conditions for which the 
appropriate boundary action reads 

Sa[^] = ±z ^ /" d'^x^/^^/^ (^_^+ + ^+^_) , (2.9) 

^ J r—ro 
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where h is the determinant of the induced metric on the boundary, and where 
the upper signs corresponds to imposing = at r = ro, while the lower 
sign corresponds to requiring 5\E'_ = 0. Which of the two Dirichlet conditions is 
imposed is a matter of convention because they can be related by inverting the 
mass M — )■ —M [24j. For definiteness, here we make the choice 

(5^+ = at r = ro. (2.10) 

Now, note that since we impose f l2.10p it is possible to add a "UV action" 
5'[/y [^^+1 on the cut-off boundary which does not obstruct the variational problem. 
If d is even, this boundary theory describes a chiral spinor, whereas for d being 
odd is a Dirac spinor. In particular, as in [T3] we shall introduce a kinetic term 
for the field and hence make it dynamical 

Suv['^+] = Z I d^x v^^+^,(r,x)^+ , (2.11) 

where Z is an arbitrary constant and l/)^{r,x) is the kinetic operator for a the- 
ory with dynamical exponent z. With F- the appropriate Dirac matrices in d 
dimensions, that is a = {t, x^, . . . ,a^_i}, this operator read^ 

lp,{r,x)=t0 = iTVd^,, (2.12) 

which becomes the usual kinetic term for a chiral fermion in the relativistic case 
z = 1. The addition of a spin-connection w^q^ on the brane is left out. In general, 
this would lead to mass terms, but for even d the spinors are chiral and the term 
with the spin-connection vanishes. For odd dimensions, this can lead to Dirac 
mass terms, but we will not consider this possibility here. 

This is how we introduce an elementary fermion in the theory. Its free dynamics 
is described by the boundary action f l2.1ip and the contribution from coupling this 
excitation to the strongly coupled CFT is described by the general relativity (GR) 
action (12. 7p and (12. 9p . evaluated on-shell. The total action 

Sfuii[^+] = Sfm + Suv[^+], (2.13) 

evaluated on shell then determines the two-point function of the elementary fermi- 
on field \Ef+, where we eliminated \Ef_ in favor of \E'+ via the Dirac equation. Note 
that we cannot add a tree level mass term on the UV brane because this would also 
require = for consistency of the variational principle. However, a boundary 
mass for the chiral fermions can be introduced through a Higgs mechanism. Finally, 

•^The sums over the indices here are, a priori, over d + 1 dimensions, including the radial 
direction. However the latter gives vanishing contribution as ^'+r-^'+ — by bulk chirality. 
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we observe that a dilaton coupling in the UV action can be included, especially in 
view of the fact that the dilaton indeed couples to the brane fields with a factor 
exp(— 0). But, the value of the dilaton at a fixed radial position r = ro can be 
absorbed into the constant Z in (12. lip . 

We also have to point out an important subtlety which has a well-known analog 
in the case of bosonic real-time correlation functions [22] . The setting above would 
correctly give rise to an effective action for \E'+ only for Euclidean signature. To 
see this subtlety, let us ignore the UV action (12. lip for the moment and set Z = 0. 
Then the prescription where we keep both terms in (12. 9 p would produce a Green's 
function for the operator that couples to the source \E'+ which is not of the correct 
form in general. In particular, it would yield a real expression that is not generally 
true for real-time correlators. The way to circumvent this problem is well-known: 
to compute the Green's function we should ignore the first term in (12. 9p . 

Namely we should replace (12. 9p by 



S9 = igf d''xV-h^^+^^ , (2.14) 

J r=ro 

where we included an additional factor of two to be consistent with the usual 
conventions. Equivalently, we can also use the prescription provided in |26| [2l]. 

Let us now become more concrete and specify the boundary space-time dimen- 
sions to be (i = 4. Also, with slight abuse of notation, we write the four- component 
Dirac spinor \Ef in terms of two-component spinors \E'+ and \E'_ as 

^ = (1:) • (2-15) 

The four- dimensional gamma matrices can be expressed in terms of a- = (1, c?) 
and a- = (—1, a) with cr- being the Pauli matrices in the following way 

r— f ° ^) ■ (2.16) 



^ 

The actions (12. lip and (I2.14p then become 



Suv['if+] = -Z I d'^a; y^^V^,(r,x)^+ , ^,(r,x) = , (2.17) 

J r=ro 

Sa = -i9f / d'^a;^^v^^-^+ , 

J r=ro 



(2.18) 

where again a = {t,Xi, . . . ,x^^i}- The two-component spinors and \E'+ are not 
independent but are related through the boundary condition at the horizon of the 
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Go. (p) = - lim r^^^^^ir, p) , (2.2I) 



black brane, which should be chosen as the in-falling one for the retarded Green's 
function. Following [26], we define Fourier-transformed spinors on each constant 
r-slice as 

/J d 
^V^±(r,p)e^P'^^^ p^ = i-oj,k). (2.19) 

The Dirac equation then imposes a relation between the chiral components tp^ 
and tp- of the form 

ip-{r,p) = -i^{r,p)'^+{r,p) , (2.20) 

where ^ is a two-by-two matrix. This matrix yields the Green's function of the 
operator 0_ in the CFT that couples to the source 4>+{p) when evaluated on the 
boundary [26] 

r—^ro 

where Go_ denotes the Fourier transform of the Green's function which takes the 
form i({C-(x), 0*_{x')})Q{t - t'). 

Next, substituting fl2.20p into the Fourier transforms of (12.171) and (12.181) . we 
find for the on-shell action (12.131) that 

5/„.[^+] = -^_ ^v^^lb/v^^(^'P) + ^^^(p)]^+' (2-22) 

where Ipzip) = ~^-^a^Pn- We regard this expression as the effective action of the 
elementary field ^+ that is coupled to an operator O- in the CFT with scaling 
dimension 

A+ = + M , (2.23) 

consistent with (II. 2p . The convention here is to write A4. for an operator (9_ 
that couples to \E'+. The effect of this coupling is described by the first term in 
(I2.22p . Next, to obtain the Green's function of a canonically normalized field, we 
transform — \&+Z'~^/^rQ^ '^''^^ and find the Green's function of the elementary 
field as 

Gr{t,,p) = - {r^oV{ro)l/>M + ^rl^'V^ro) e(ro,p))'' , (2.24) 

where we remind that the blackness function V{r) appears in the metric given by 
(12. 2p . In the limit of zero temperature, or for temperatures small compared to 
the UV cut-off, we have rQ/r^ ^ I so that we can approximate V^tq) ~ 1. Once 
(I2.24P is obtained in this way, we take a double scaling limit 

ro ^00, gf^ 0, gfrl'^^~^^'^ = const. (2.25) 
Furthermore, we note the following consistency checks: 
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(i) In the case of g/ = and finite tq, hence the dynamical source field is 
decoupled from the CFT, we find the correct propagator for the free field 



(ii) In the case of non-vanishing Z we calculate the correlation function by tak- 
ing functional derivatives of the generating function Z[J^] where J_ is the 
Legendre transform of \I'_|_. Schematically, this reads 



G^oc^^ with Z[J_] = / Pvl/+e^^/-[^+l+^/^'^-(H-^-+^-*+), (2.26) 



where Sfuu is given by f l2.22p . This provides an alternative way to arrive at 
the result fl2.24p . At the end of this calculation one takes the limit fl2.25p as 
explained above. 

(iii) In the case of Z = 0, that is the UV dynamics for the source \E'+ is turned 



inverse of the result in the standard quantization. This is expected because 
in the alternative quantization above, we integrate over the source, see [27], 
rather than fixing it. 

In passing, let us comment on the range of the parameter M for general z. 
In analogy with the AdS case, 2; = 1, we find that, among the possible bound- 
ary asymptotics of \E', \E'+ is non-normalizable near the boundary for M > z/2. 
Therefore it should correspond to the source term in this range. On the other 
hand, in the range < M < z/2 both \Ef_|_ and are normalizable, therefore 
one has the possibility of choosing either of them as the source. This corresponds 
to two different quantizations [27J. Using the fact that the system is invariant 
under the exchange \E'± — )■ and M — )■ — M, one instead extends the range of 
M to M > —z/2 and always fixes \E'+ to be the source. In this way, one covers all 
possible quantizations. In particular, the range —z/2<M<0 now corresponds 
to the alternative quantization where \E'_ is chosen as the source. This is what we 
shall do in the following, i.e., we will always choose \E'_(. as the source but consider 
the wider range M > —z/2. 

As in the AdS case [24], we can easily derive a first-order differential equation 
that determines the matrix ^. For this purpose, let us specify to the case of d = 4 
for which \E^_|_ describes a two-component spinor. A similar analysis can be done 
for d = 3, as long as z < d — 1. The up and down components u± and d± of \E'± 
are defined as follows 



+ 





(2.27) 



'*We thank Hong Liu for a discussion on this issue. 
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where the overall r-dependent factor is introduced to simplify the fluctuation equa- 
tions. Then, the eigenvalues of the matrix ^ above are given by the ratios 



^+ = t— , ^- = ':r ■ 2.28 

Employing the rotational symmetry of the problem, we can choose a basis where 
the particle moves along the z-axis, i.e. k = (0, 0, k). Then we first derive the fol- 
lowing first-order decoupled equations, that relate the plus and minus components 
in (K27} as 



i{oj + k)u+ = A{~M)u_ , i{uj - k)u^ = A{M)u^ 

i{u - k)d+ = A{-M)d^ , i{u + k)d^ = A{M)d^ 



(2.29) 



where we introduced the shorthand notations 



A{M) = r{rVdr - M) . (2.30) 



Using the first-order equations above we can now derive a first-order equation for 
^+ and ^_ defined in fl2.28p as 

r'^VdrU + 2Mr^± = -cl; ^ A; + (-w ± k)il . (2.31) 

The in-falling boundary conditions at the horizon correspond to 

i±{rh,p)=i. (2.32) 

Once the function ,^-|- is obtained, the Green's function for the elementary field 
coupled to the CFT is given by the equation (12.241) above. In passing, let us also 
observe that by a change of variables we find that ^±{r,uj, k) can be expressed in 
terms of the following ratios 

=«4pf)=«H^'^) . (2.33) 

3 Results 

3.1 The relativistic CFT with z = \ 

In the case that the elementary fermion is coupled to an operator 0_ in a 
relativistic CFT, the background is given by an AdS black brane, the case z = 1 in 
(12. 2p . For simplicity of the notation, we choose again the two-component spinors 
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introduced in fl2.27p . The kinetic term f l2.12p in the zero-temperature hmitH and 
for z = 1, then takes the form 

, — * — * 

rQiPi{rQ,u, k) = u — a ■ k . (3.1) 

Let us now focus on the case A; = 0. The solution to the differential equation fl2.3ip 
with infalling boundary conditions ^±(r = 0) = +i reads 

J,, ■ {i^)+e"'"+ilj 1^) 

where Ja{x) are Bessel functions of the first kind. Focussing on the case — 1/2 < 
M < +1/2, performing an expansion for large r and keeping only the leading order 
term, we find the following expression [21] 

e(ro, c.) ^ - (2ro)-^^^ l^^'+M) ^'^^^^^^^^^"^ • (3-3) 

The bulk mass M and the scale dimension of the operator O- is given by A+ = 
d/2 + M. From f l2.24p we then find the Green's function for the up (-I-) and down 
(-) components of the spinor which reads 

Gr{6,u) = — , (3.4) 

where we remind the reader that in (13.41) a two-by-two identity matrix is under- 
stood, and where we defined the constant 

Note that in case qm vanishes, hence is decoupled from the CFT, one should 
add a factor of ie with e > in the denominator of (13. 4p . as usual for the retarded 
Green's function. Furthermore, qm is positive definite in the range —1/2 < M < 
+1/2 provided Qf/Z is positive definite. 

As can be seen from (13. 4p . the self-energy is proportional to u'^'^'^ which for 
M < diverges in the infrared. This is an interesting situation with possible 
applications in condensed-matter physics. Note also that the retarded Green's 
function satisfies 

G^^i6,u) = -GRi6,-u), (3.6) 



^As noted in the introduction, we obtain the zero T resuh, by first performing the calculation 
at infinitesimally small horizon and then take the limit T — > 0. 



14 



which, using fl2.3ip . can be generahzes to non- vanishing spatial momentum as 



Tr w) = -Tr Gr(A?, -tu) . (3.7) 

Now, from the expression (13 ■4p we find a pole at a; = 0, and a branch-cut 
that we can place on the negative imaginary axisjf] The location of the only other 
possible pole in this expression depends on the sign of g. More concretely, in the 
range —1/2 < M < 1/2 there is no pole in the upper half-plane of the principal 
sheet provided that 

> . (3.8) 

Thus, the expression (13. 4p satisfies the Kramers-Kronig relations. Moreover, using 
Cauchy's theorem for a contour in the upper half plane with a semi circle closing 
at infinity, we find 

/oo 
duTTGR{k,u) = 2TTi . (3.9) 
'OO 

For this identity, it is essential that M < 1/2 such that at infinity the fall-off of 
the Green's function goes like u~^, and such that the pole at w = does not 
contribute. A sum rule now follows straightforwardly. Introducing the spectral 
density functionlll 

p{k, cu) = ^lm Tr [Gnik, u)] , (3.10) 
and employing (13. 6p as well as (13. 9p . it then follows thatlf] 

/oo 
dup{k,u) = l. (3.11) 
-oo 

This is the ARPES sum rule required for one-particle fermionic states. It holds 
for any value —1/2 < M < 1/2. The fact that the sum-rule is violated outside 
this range is in accord with the observation that M = —1/2 corresponds to the 
unitarity bound for fermions in a relativistic CFT. 

For non-vanishing momenta, one can easily redo the calculation. The non- 
trivial part of (I2.24P is the self-energy that is given in terms of C,{r,p). This is 
of course equivalent to knowing the two-point function of the operator O- in the 

^There is a physical reason for this. As one turns on temperature the branch-cut disintegrates 
into an infinite series of poles which should appear on the negative imaginary axis [5] . We thank 
David Vegh for a discussion on this point. 

^Note the relative minus sign in this definition with respect to the more conventional definition, 
which we denoted by G in the introduction. This is because we define our correlators with an 
overall minus sign, following the convention of 1241 . 

^Employing p.4p for the Green's function in p.lOp . we also checked analytically that the sum 
rule is obeyed. 
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CFT, which is a known result in the hterature also for non-zero spatial momentum 
[3T| l32| l33t [21] . In the two-component notation, the result can be expressed as 



Go_{p)=T{p)a^5/p^, (3.12) 
where the function T(p) for large values of r is given by 

T{p) ~ (2r) 



-2M n\-M) 

va + M) 



^p2M-lg-m(Af+|) ^ ^^2_ |^|2 ^ ^ > +1^1 ^ 

^p2M-lg+i7r(Af+i) ^ p ^ a/u;2 _ |^|2 ^ W < , 



X < 

In order to investigate the analytic structure, we have to find the analytic contin- 
uation of the function T{p) in the entire complex uj plane. We note that to obtain 
the second line in (13.131) with u < —\k\ from the first with uj > +\k\, we take 
p — )■ e^'^p. Similarly, to obtain the last line we take p — )■ e^'^/'^p. Therefore, the 
single expression 

T{p) = (2r)~^^ |;|^^| e-(^+iy^-^ , p ^ - W , (3-14) 

covers the entire complex u plane provided we include a branch-cut that runs 
from a; = to a; = —ioo and take the first sheet to make it single-valued. So in 
particular, this means that —1 will be represented by e~^^'". Now, the full Green's 
function (12.241) is given by 

Gn{k, u) = ^ (u + a-k) , (3.15) 

p2 (l~gMe-'''(^^+-2^p^M-i 



where the constant qm is defined in (13.51) . The poles of the retarded Green's 
function coming from the first term in the denominator of (I3.15P are at 

uj = ±k, (3.16) 

and there are no poles in the upper half-plane of the principal sheet if 

^Af > . (3.17) 

Thus, the Kramers- Kronig relations are obeyed. Also, using the symmetry prop- 
erty (13. 7p and the fact that the large frequency behavior is the same as for k = 0, 
one finds again that (13. 9p is satisfied, and hence the sum rule is obeyedj^ 



^Using (13.121) and (I3.13P in the expression for the full Green's function in (|3.10p . we also 
checked analytically that the sum rule for non- vanishing k is obeyed. 



16 



3.2 The Lifshitz case 



For values of the dynamical exponent z different from one, there are a number of 
important changes. We consider again the case of vanishing temperature, implying 
V^(r) = 1, and first assume both k and w to be non-zero. The kinetic term in fl2.22p 
then reads 

rlJj)X^,k)=u}-rl-^o-k . (3.18) 
After rescaling ^+ — )■ \E'+Z^^/^rQ^ the action fl2.22p becomes 

^ -^^i^i[%rl^'ar,p)+u~rr'a-k]i^+. (3.19) 

Now, at the end of our calculation, we will remove the UV cut-off by taking the 
limit ro — )■ oo. Consequently, a possible divergence in the first term in fl3.19p should 
be absorbed in the redefinition of the parameter Qf. However, the divergence in 
the last term should be cancelled by a counter-term action, which is parallel to 
what happens in field theory. To explain that point, let us consider a single spin- 
component ip of a. fermion in a Lifshitz invariant free theory. The kinetic term 
reads Skin oc J du d'^~^k il)*{uj + rik^)ilj. The coupling rj is classically marginal, 
with the scaling u oc k^ and the classical scaling dimension for ip being {d — l)/2. 
In presence of a UV cut-off in spatial momentum of order oc ro, however, 
the Lifshitz scaling is broken and one has to consider the contribution of relevant 
terms. The kinetic action then becomes Skin oc J du d'^^^k iIj*{u} + rjk^ + 'fjk)ip 
where fj scales like fj oc A^~^, hence the term is relevant in the IR limit. If we want 
to maintain Lifshitz scaling in the IR, we have to renormalize and remove fjk, and 
we are left only with a classically marginal spatial term rjk^. In the end of this 
procedure, from f l3.19p we obtain the following Green's function 

GR{k,u) = -(u + va-k k^-' + I r^,+'ar, p))'' ■ (3.20) 

We also note that the matrix C,{rQ,p) in f l3.20p is not arbitrary but its form is 
determined by the scale dimension of the operator 0_ and the Lifshitz scaling. 
Choosing again a basis where the particle moves along the 2;-axis, we can write 
k = (0, 0, k) and obtain 

.i,„^r-«,.„.p)^.-|,(^)^.^|,(^). (3,21) 



^■0 



Here we indicated explicitly that the matrices C,i 2 are only a function of the ratio 
u/k^. This is of course expected as it is essentially the Green's function of the 
operator O- in the Lifshitz CFT. We can see this in the holographic picture directly 
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from equation (12.311) . First consider T — )■ which means setting V{r) = 1. Now 
we perform a change of variables x = r/k under which (I2.3ip becomes 



dx 



y^{x) + 2Mxi^{x) = ^ T 1 + ± l) , (3.22) 



where u = u/k^. Note that uj and k only appear in the combination cJ, and the 
boundary condition (I2.32p also respects this. Then, the scaling in (13.211) follows. 

Unfortunately, for arbitrary z we cannot solve (12.311) analytically when both uj 
and k are nonzero (an exception is the special case M = and z = 2, see [3^). 
However, for the particular case = and a; 7^ the solution to (I2.3ip for infalling 
boundary conditions (I2.32p reads 

(^)+e-(^+i)j_M^i (^) 

«---'- ,:::(i)-e-.f^^>;::(i) - 

— +2 Vzr^; - — ~2 \zr^' 

where Ja{x) are again Bessel functions of the first kind and where a two-by-two 
identity matrix is again understood. Its expansion for large values of r can be 
obtained as 

P / 1 M\ 

e(r,u;) -r-^' \)l ^ e-M.+^)(2^)--a;- . (3.24) 

Note that these expressions hold for generic values of M, however for the special 
values M G Z + I logarithmic terms appear. Furthermore, in the situation of 
a; = and A; 7^ the solution to (I2.3ip reads as follows 



^+Af 1 (--) 
A+A^+i (--) 



(3.25) 



where /^(x) and Ka{x) are the modified Bessel functions of the first and second 
kind, respectively. Also, the infalling boundary conditions have to be modified to 
^±(0, k) = +1, and the expansion for large values of r in both cases reads 

e(r, k) ^ -(2r)-^^^|^ P^-' k-a. (3.26) 

Again, these results hold for generic M, but for M G z(Z + |) logarithmic terms 
appear. 
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In order to investigate the analytic structure of the Green's function f l3.20p . 
let us consider the general case u ^ Q ^ k. Since we do not know the analytic 
expression for ^{r,p) in this case, we have to make the assumption that all possible 
non-analyticity of f l3.20p is in the lower half-plane, i.e., Kramers-Kronig relations 
are satisfiedo With this assumption, we can evaluate the integral 

dujp{k,uj) = — dujTi[GR{k,uj)] , (3.27) 

where we employed the definition of the spectral density function given in f l3.10p . 
In particular, we can deform the integration contour from the real axis to the 
infinite semi-circle in the upper half plane, i.e., we write u = Re^^ with < < vr 
and send i? — > oo. Let us consider the second expression for the matrix ^{r,p) in 
fl3.2ip . The value of ^2 on this semi-circle can easily be determined by noting that 
the limit u ^ 00 with k finite is the same as — )■ and u finite. Using f l3.24p this 
gives lim^_j,oo ^2 — const, on the infinite semi-circle. Let us denote this semi-circle 
as C. Then, we find 

/ dujp{k,uj) = — [ dujTT\GR{k,uj)] =1 , (3.28) 
i-00 2z7r Jc 

iff M < z/2, regardless the value of k (as long as it is finite). In this case, also the 
pole at a; = does not contribute. Therefore we again find that the sum-rule is 
obeyed for —z/2 < M < z/2. This generalizes the result in the previous section 
for z = 1. 

In this general case we simply assumed that the Kramers-Kronig relation holds 
and then showed that the sum-rule is obeyed in the allowed range of M. For the 
particular case of = 0, however, we can also the Kramers-Kronig analytically. 
This is a straightforward generalization of discussion below equation (13. 4p where 
for generic z, in place of (13.41) we have 

Gr{0, uj) = : L , ,M^i, • (3.29) 

The constant (13. 5p in the present case generalizes as 



= f (2.)-^ ir)f7^r-««+' , (3.30) 



which is positive definite in the range —z/2 < M < z/2 for Qf/Z being positive 
definite. We again find a pole at w = 0, a branch-cut on the lower imaginary axis 
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We study the analytic structure for this general case and at finite temperature numerically 
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and there is no pole in the upper half-plane of the principal sheet provided that 

9M,z > . (3.31) 

Clearly, all non-analytic behavior is in the lower half-plane and the Kramers- 
Kronig relations are satisfied. For negative mass values, the self-energy is infrared 
divergent, which signals an IR singularity that may have interesting applications 
for condensed matter physics. 

One can also ask what happens beyond the range \M\ < z/2. Explicit calcu- 
lations show that for M > z/2 the requirement of analyticity in the upper 
half plane can only be satisfied when g/fZ < 0. However the sum-rule integral 
fl3.28p then yields —1, which signals a violation of unitarity and which is indeed 
consistent with the "wrong" choice for the sign of gj/Z. On the other hand, in 
the opposite range M < —z/2 one finds that there always is a singularity in the 
upper half plane, as a result of which both unitarity and causality is violated. We 
leave the details of this calculation to [391. 



4 Discussion 

In this work we introduced a method to compute correlation functions of elemen- 
tary fermion fields that satisfy canonical commutation relations, within a strongly 
interacting CFT with an arbitrary dynamical scaling exponent z, in the context of 
the holographic correspondence. Our emphasis is to satisfy the sum-rules that are 
obeyed by these elementary fermion fields, as observed in ARPES experiments. 
We determined the condition for satisfying the ARPES sum-rules imposed on the 
mass of the dual bulk fermion field as —z/2 < M < z/2 for an arbitrary dynamical 
exponent z. In the field theory this requirement corresponds to coupling the ele- 
mentary field to an operator (9„ with scaling dimension above the unitarity bound 
and relevant in the IR. 

Let us expand a little more on the latter condition. It is clear from (13.291) that 
demanding the interaction term be relevant in the IR is M < z/2. One can also 
check this directly in the corresponding field theory action: Such an interaction is 
given by J dt cfi~^x g \E'_|_(9_ in the action. The question is whether the coupling g is 
irrelevant or not in the UV. The scale transformation is x — hr^x, and t — )■ hr^t. 
The weights of the elementary fermion \E'+, the operator and the volume term 
under this transformation are {d — l)/2, A+ and —{d — 1 + z), respectively. Then 
the dimensionless coupling that one constructs from g is g = gA^+~^ ~. For the 
beta-function of this coupling to be negative, i.e. being relevant in the IR, one 
should have A+ < z + which indeed corresponds to M < 2;/2 using (II. 2p . 
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As a by-product of our analysis, we present analytic expressions for the fermion 
correlation functions in strongly interacting non- relativistic CFTs@ 

There are various directions one can extend our work. The most immediate 
generalization involves turning on nonzero temperature and chemical potential 
P^ . The zeroth sum rule fl3.9p should also be satisfied in the case of nonzero 
temperature. This will be a non-trivial test of the generalization of our method 
to arbitrary temperatures. This step will be crucial of one desires to compare the 
results of holography with real ARPES data for systems with similar properties 
such as bi-layer graphene. 

Another question involves coupling the elementary field to more than one CFT 
operator. Indeed, this should be a more general situation in the ARPES experi- 
ments: the excited elementary fermion may interact with the strongly coupled CFT 
through many operators. This case was considered in [16] in a semi-holographic 
fashion. In order to study this situation in a more "holographic" manner, one can 
imagine turning on a source for each CFT operator O^. In addition, if one 
turns on another elementary fermion, say X- the UV cut-off surface, couple it 
to all of the sources on the boundary, and calculate the two-point function of X- 
instead of ip^, one may be able to achieve this goal. 
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